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We theoretically investigate the response of optical phase 
conjugators to incident probe pulses. In the stable (sub- 
threshold) operating regime of an optical phase conjugator 
it is possible to transmit probe pulses with a superluminally 
advanced peak, whereas conjugate reflection is always sublu- 
minal. In the unstable (above-threshold) regime, superlumi- 
nal response occurs both in reflection and in transmission, at 
times preceding the onset of exponential growth due to the 
instability. 
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I. INTRODUCTION 

A variety of mechanisms is now known to give rise 
to superluminal (faster-than-c) group velocities, which 
express the peak advancement of electromagnetic pulses 
reshaped by material media: 

1. Near-resonant absorption anomalous disper- 
sion in the linear regime of an absorbing medium 
forms the basis for this superluminal reshaping 
mechanism. 

2. Reduced transmission or evanescent wave forma- 
tion (tunneling) in passive dielectric structures 
[|| ||: this reshaping mechanism has been at- 
tributed to interference between multiply-reflected 
propagating pulse components in the structure || . 

3. Soliton propagation in dissipative nonlinear struc- 
tures 0: superluminal group velocities can oc- 
cur in such systems via nonlinear three-wave ex- 
changes, as in stimulated Brillouin backscattering 
in the presence of dissipation. They also occur in a 
nonlinear laser amplifier ||. 

4. Pulse propagation in transparent (non-resonant) 
amplifying media [pyTcfl. Superluminal pulse re- 
shaping in this regime has been attributed to either 
the dispersion || or the boundary reflections |l(]] 
of the amplifying medium. 

5. Tachyonic dispersion in inverted two-level media: 
the dispersion in such inverted collective systems 
is analogous to the tachyonic dispersion exhibited 
by a Klein-Gordon particle with imaginary mass 
[pdf . Consequently, it has been suggested Jl^] 



that probe pulses in such media can exhibit su- 
perluminal group velocities provided they are spec- 
trally narrow. Gain and loss have been assumed to 
be detrimental for such reshaping. We note that 
Ref. p2| describes an infinite medium and bound- 
ary effects on the reshaping have not been consid- 
ered. 

In this paper we establish a connection between op- 
tical phase conjugation and tachyonic-type behavior in 
a finite medium. To that end we study the reflection 
and transmission of pulses at a phase-conjugating mirror 
(PCM), both in its stable and in its unstable operat- 
ing regimes. Similar features in other parametric pro- 
cesses, such as stimulated Raman scattering and para- 
metric down-conversion will be treated elsewhere. The 
dispersion relation in infinite PCM media, which allows 
for superluminal group velocities, was derived by Lenstra 
plU . Here we address the questions: can superluminal 
features be observed in the response of a finite PCM to 
an incident probe pulse, and how can they be reconciled 
with causality? 

Pulse reshaping by a PCM has been studied before 
JLlj , but never in the context of superluminal behavior. 
We show (Sec. OB) that in the well-known stable oper- 
ating regime of a PCM, superluminal peak advancement 
can occur in transmission but not in conjugate reflection, 
as was briefl y pr esented in Ref. In the unstable 

regime (Sec. [I C), we demonstrate the existence of su- 
perluminal peak advancement in both the reflected and 
transmitted wave response to spectrally narrow analytic 
probes. Further insight into these processes is gained by 
spatial wavepacket analysis (Sec. III). For sharply-timed 
(abruptly modulated, non-analytic) signals (Sec. IV) the 
occurrence of spectral components in the frequency zone 
where the gain is highest cannot be avoided. These com- 
ponents trigger early exponential growth, which is inde- 
pendent of the probe pulse shape. In all cases we develop 
criteria for the optimization of probe pulse shapes, aimed 
at maximizing their superluminal peak advancement and, 
in the unstable regime, the delay of the instability onset. 

Our phase-conjugating mirror consists of a nonlinear 
optical medium with a large third-order susceptibility 
X^ 3 \ confined to a cell of length L. Phase conjugation 
is obtained through a four-wave mixing (FWM) process 
p6[ : when the medium is pumped by two intense coun- 
terpropagating laser beams of frequency too and a weak 
probe beam of frequency lu p = ujq + <5 is incident on the 
cell, a fourth beam will be generated due to the non- 
linear polarization of the medium. This conjugate wave 
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propagates with frequency uj c — loq — 6 in the opposite 
direction as the probe beam (see Fig. |l|). 
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FIG. 1. Phase conjugation by four-wave mixing. Solid 
(dashed, dotted) arrows denote pump (probe, conjugate) 
beams. 



II. PULSE REFLECTION AND TRANSMISSION 
TEMPORAL ANALYSIS 



A. Basic Analysis 



The basic semiclassical one- 
scribing this FWM process 
tuting the total field E(x, t) 

refer to the two pump beams 
conjugate beams respectively) 
a nonmagnetic, nondispersive 
a nonlinear polarization 



-dimensional equations de- 
are obtained by substi- 

= E^iAp/^.') = 
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Selecting the phase-conjugation terms for P/vi(a;,i) = 
X^ 3 'E 3 (x,t), assuming the pump beams to be non- 
depleted and applying the slowly varying envelope ap- 
proximation (SVEA) results in |lq| 
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Here f^*^(a;,i) denotes the complex amplitude of the 



probe (conjugate) field and k 



is the coupling strength (per unit length) between the 
probe and conjugate wave. The dispersion relation for 
an electromagnetic excitation in this pumped nonlinear 
medium is given by [|13[ 
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where S 
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and k 2 is the squared 



The group velocity 35/ dk in the medium can be shown 
from (|^) to be always larger than c, the speed of light in 
vacuum, and the dispersion relation is therefore of the 
tachyonic type (see Fig. §(a)), analogous to the one for 
inverted atoms |12| . Because of the superluminal group 
velocity caused by this dispersion, the question arises 
how wave packets will be reshaped by a PCM. In the 
early 80's, Fisher et al. studied the phase-conjugate re- 
flection of pulses of arbitrary shape at a PCM No 
tachyonic effects were found. They distinguished between 
two different operating regimes of the mirror, stable (if 
k L < |) and unstable (if k L > §)■ In the stable 
regime, the response of the mirror is always finite and 
scales with the probe input. The unstable regime cor- 
responds to self-generation of conjugate reflection from 
arbitrarily small probe input, followed by exponential 
growth until saturation is reached (due to depletion) . Re- 
cently we predicted the occurrence of superluminal ad- 
vancement of the peak of a suitably chosen input pulse 
upon transmission through a stable PCM p5| . Here we 
present the full analysis of both the stable and the un- 
stable regimes and demonstrate the existence of superlu- 
minal effects in both. 

In order to study pulse reflection and transmission at 
a PCM, we use the two-sided Laplace transform (TSLT) 
technique introduced by Fisher et al. [llif . The basic 
approach is summarized in App. A. To begin with, the 
reflection and transmission amplitudes for a monochro- 
matic pro be be am incident on a PCM are given by 
(Eqs. ([A2]) and QAjj), at x — and x — L, respectively) 
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Now consider a probe field £ p (0,t) incident on the PCM 
at x = 0, which satisfies the basic TSLT premise that it 
decreases faster than exponentially as t — > — oo We 
obtain the expressions for the resulting reflected phase- 
conjugate pulse at x = 0, and the transmitted probe 
pulse at x = L from the inverse TSLT (see App. A). 
In taking the inverse TSLT we separate the singularities 
of r cp {is) and t pp (is) in the right-half s-plane, which give 
rise to unstable exponentially growing solutions, from the 
singularities in the left-half s-plane, which correspond to 
stable solutions. The final expressions for < n^L < 
are 



£c(0,t) 



wavevector of both the probe and the conjugate waves. 
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and X\ the nontrivial solution of sin(X) = ±X/(kq L). 
In the second terms in ([?]) and (||) si is the unstable pole, 
and hi and fa are the residues of the reflection and trans- 
mission amplitudes at this pole. After rewriting the first 
term in ([?]) and (|8|) as a Fourier transform, it is straight- 
forward to analyze £ c (0,t) and £ p (L,t) numerically. If 
KoL < 7r/2 there are no singularities with Re(s) > 0, so 
hi = fa = and the second term in ()?]) and (||) does 
not contribute. Hence, the stable regime is defined by 
kqL < ir/2, and the unstable by kqL > ir/2. We now 
start by studying pulse reshaping by a stable PCM, and 
then move on to the unstable regime. 



B. Stable regime 

Fisher et al. jlj] analyzed the phase-conjugate reflec- 
tion at a PCM in the stable regime (koL < ir/2). For a 
gaussian input pulse they found a delay of the peak of 
£ c (0,t) with respect to that of £ p (0,t). Here we wish to 
emphasize that a suitably chosen incident pulse in this 
stable regime is reshaped in such a way that its peak 
emerges from the cell before the time it takes to travel the 
same distance in vacuum. The central result, obtained 
from Eqs. (^) and (||) with h\ — h 2 — 0, is depicted in 
Fig. |(b). 
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FIG. 2. (a) Dispersion relation (solid line) in a PCM. 
The dashed (dotted) lines correspond to the dispersion re- 
lation for the probe (p) and conjugate (c) waves in vac- 
uum and fco = i^o/c. The marked positions 1-3 in- 
dicate the frequency components in the incident pulses 
which give rise to the transmitted curves 1-3 in (b). (b) 
Transmitted probe pulse \£ p (L,t)\ at x = L for an in- 
coming gaussian £ p (0, t) — e~ at e lSat at x = (thick 
solid line). The temporal width (FWHM) of the incom- 
ing pulse A t = 2-^/ln 2/q = 3.3 L/c, the spectral width 
A,5 = 4V q In 2 = 1.7 c/L and kqL = 1.4. The vertical line 
indicates the time tt r ~ L/c needed to traverse the cell in 
vacuum. 

The condition on the frequencies in the input signal 
for observing the superluminal effect is that they should 
be centered away from the gap in the dispersion rela- 
tion. The reason for this is that a pulse which is cen- 
tered around Sq — contains positive as well as neg- 
ative frequency components whose respective positive 
and negative group velocities interact and compensate 
in such a way that no superluminal advancement occurs. 
For a pulse centered further up on the dispersion curve 
(Fig. |2|(a)), however, superluminal peak advancement is 
obtained (at S > 2 c/L, curves 2 and 3 in Fig. ||(b)). 
This superluminal peak advancement depends on three 
parameters: the temporal width of the incoming pulse 
A t , the central frequency of the incoming pulse 5q and 
the coupling strength in the medium kq. By fixing kq 
and varying simultaneously A t and So the superluminal 



3 



effect can be (numericalfy) optimized in several ways. 
In absolute terms, we find a maximal attainable peak 
advancement of ~ 0.88 L/c. But since this advance- 
ment, even though large, would only be a small effect if 
the pulse is broad in time, it is also useful to optimize 
the ratio r — peak advancement/pulse width. We ob- 
tain a maximal relative peak advancement of r ~ 0.08. 
In phase-conjugate reflection no such superluminal effect 
appears: the time at which the peak of the reflected sig- 
nal emerges at x = is always later than t = 0, the time 
at which the maximum of the input pulse entered the 
cell. 



C. Unstable regime 

We now move on to the unstable regime (kqL > tt/2), 
for which Eqs. (0) and (g) have /ii,/i2 7^ 0. Figure ^ 
shows the reflected phase-conjugate pulse in this regime 
for an incident gaussian centered around frequency Sq. 




FIG. 3. Phase-conjugate reflected pulse |£ c (0, t)\ at 
a PCM in its unstable operating regime (koL — 1.7) as 
a function of time (in units of L/c) for incident gaus- 
sian-shaped pulses £ p (0,t) = e~ at e' Sot with different So (in 
units of c/L). Curves a, b and c correspond, respectively, to 
S = 0.1, 0.28 and 0.6. |£p(0,t)| (thick solid line) has tem- 
poral width A t = 2^/ln2/Q = 24 L/c and spectral width 
A s = 4v / oTn2 = 0.23 c/L. 

For #o = 0.1 c/L we see that the reflected signal starts 
growing exponentially as soon as the incoming pulse 
reaches the cell. However, for an incident pulse centered 
around a frequency further away from the gap in the dis- 
persion relation (at <5 ~ 0.28 c/L), the reflected pulse 
exhibits a local maximum before the exponential growth 
sets in. This peak is clearly advanced with respect to the 
peak of the incoming signal. Since KqL is chosen close to 
7r/2, the reflected pulse is greatly amplified pl| . 

For large So, where the dispersion relation becomes 



asymptotically linear (8 — ±c(fco ± k)), the reshaping 
of the reflected pulse is only minor and the superluminal 
effect is no longer noticeable for this pulse. In order to 
optimize the superluminal response preceding exponen- 
tial growth, one must thus have So small (for maximum 
advancement) and KoL close to tt/2 (for maximum inten- 
sity and delay of exponential growth). 

Just as in the stable regime, we need to optimize three 
parameters simultaneously: kq, 5q and A t . One way of 
doing this is by using the exact expression (Q) and ana- 
lyzing it numerically, but this does not give much insight 
into the interplay between these parameters. Another 
method is to find an approximation of the exact result, 
which allows for an analytical treatment of a certain class 
of incident pulses and yields quantitatively good agree- 
ment with the exact result in that case. 

In order to obtain such an approximation we con- 
sider the conjugate reflection amplitude (Q) in the limit 
S -C k c, corresponding to incident pulses with a large 
temporal bandwidth compared to L / c. Equation (J4j) then 
reduces to 
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The reflected pulse becomes in this approximation 
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£ c (0,t) 



dt£* p {Q 1 t)e~ (16) 



and we see that the growth in the unstable regime t m < 

t 

behaves as ei*™i . For spectrally narrow pulses and large 
Ko, the difference between ([l6]) and the exact numerical 
result is found to be < 10%. 

Taking the derivative of |£ c (0,£)| from Eq. ( [l6| ) with 
respect to t and equating it to zero gives a condition on 
the times at which the reflected pulse is maximal (t max ) 
or minimal (t min ), 

X 2 (t) + Y 2 (t) - t m [X{t)x\t) + Y{t)Y\t)] = 0, (17) 



with 



X(t) 



dt'Re(£*(0,t')) 



,t lu 



Y(t)= / dtlm(£;(0,t))e 



t /t„ 



(18) 
(19) 



The optimal superluminal effect is found if the reflected 
intensity is large at t max and close to zero at t m i n and the 
separation f TO j„ — t max is as large as possible. In order to 
obtain the maximal relative advancement we numerically 
scan through the three-parameter space (no, Sq, At), fix- 
ing KqL and varying the other two parameters simultane- 
ously. Using (fijl), this reveals that the reflected pulse is 
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very sensitive to 5q (for fixed kq and A t , see also Fig. ||) 
and only arises for incident pulses that are sufficiently 
broad in time (A t > 20 L/c). For signals with broad 
spectra, the strong influence of frequency components in 
the instability gap prevents the formation of a discernible 
pulse response before exponential growth sets in. Fur- 
thermore, KoL should be close to tt/2, because otherwise 
the fast onset of exponential growth masks the reflected 
pulse. In the optimal case one can find an advancement 
of the peak of t max ~ 10 L/c for a pulse of temporal 
width A t ~ 25 L/c. The peak intensity \£ c (0, t max )\ ~ 5, 
\£ c (0,t min ) \ ~ 0.1 and t min - t max ~ 6 L/c. 

The results for transmission of a gaussian through an 
active PCM are qualitatively the same as for phase- 
conjugate reflection. The amplitude of the advanced 
transmitted response is different, but not the values of 
(5o and A t for which it arises. The approximation (|l4| ) 
cannot be used to describe the transmission in the sta- 
ble regime, which occurs for larger values of Sq than the 
superluminal response in the unstable regime. The as- 
sumption S <C kq c is not valid in that case. 



III. PULSE REFLECTION AND TRANSMISSION 
- SPATIAL ANALYSIS 

To gain further insight into pulse reshaping by a PCM, 
we follow an incoming gaussian pulse in space. One can 
then observe the following stages: (1) the probe wave 
packet approaches the cell; (2) it propagates as an "op- 
tical quasiparticle" (consisting of a probelike part travel- 
ing to the right and a conjugatelike part traveling to the 
left) jl9| in the cell and (3) the reflected phase-conjugate 
and transmitted probe packets leave the cell. We employ 
again the TSLT of Sec. [J. The reflected phase-conjugate 
and transmitted probe pulses at position x in the nonlin- 
ear medium and time t are given by (App. A) 
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where h r (x,S) and h t (x,5) are given by Eqs. (A2) and 
(A3). For the setup of Fig. [l], the total probe pulse 
£ p (x, t) now consists of an incoming probe pulse in the 
region x < 0, the probelike part of an "optical quasipar- 
ticle" in the cell (0 < x < L) and a transmitted probe 
pulse for x > L. The result is 



£ p (x,t)=6{-x) £^ M (0,i--) + 



PCM, 



(x-L)£^(L,t 
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(e(x)-e(x-L))£^ CM (x,t) + 
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PCM (M+ - ) + 

c 

{6{x) - 6{x - L)) £^ CM (x,t). 
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We subtract from (|22j) and (|23|) again the contribution of 
the singularities which give rise to exponentially growing 
solutions, as in Eqs. (0) and (||), and rewrite them as 
Fourier integrals over frequency S. We then consider an 
incoming gaussian and analyze |£ p (o;,i)| 2 and |£ c (a:,£)| 2 
as a function of x (and t) numerically, especially focusing 
on the intruiging pulse reshaping effects found before: the 
possibility for superluminal peak traversal times in probe 
transmission in the stable regime and the superluminal 
pulse response in the unstably operating PCM. 

Figures |I] and || show the time evolution of a probe 
wave packet incident upon a PCM. The incoming pulse 
\£ p (x, t)\ 2 is centered around frequency Sq and its width is 
given by (FWHM) As = 2 \j2a In 2, which corresponds to 
a spatial width A x — c ^/21n2/a. For clarity — \£ c (x, t)\ 2 
is plotted along the vertical axis instead of +\£ c (x,t)\ 2 . 
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where is the Heaviside step function, and similarly, 



FIG. 4. Stable regime (below-threshold) reflection and 
transmission of an incoming probe pulse at a PCM of length 
L (indicated by shaded area) . The solid (dashed) curve in the 
upper (lower) half of each time plot shows the probe (conju- 
gate) pulse |£ p (a;,i)| 2 (- \£ c (x,t)\ 2 ), which is moving to the 
right (left). Parameters used are So = 3 c/L, A x — 2.1 L, 
As — 1.3 c/L, KoL — 1.4. The dashed vertical line for 
t = 5 L/c indicates the position where the peak of the pulse 
would have been in the absence of the PCM. 

Figure || depicts the probe transmission and conjugate 
reflection in the stable regime (kqL < 7r/2) for an in- 
coming pulse of spatial width 2.1 L, centered around fre- 
quency 5q = 3 c/L in the frequency domain. These pa- 
rameters are comparable to those for curve 2 in Fig. @, 
but since we now consider \£ p (x,t)\ 2 instead of \£ p (x,t)\, 
the normalization is different. We see how the probe 
pulse approaches the cell at t — — 2 L/c. At t = 0, 
when the forward tail of the pulse has entered the cell, a 
small reflected phase-conjugate pulse has developed and 
is traveling simultaneously in the opposite direction. At 
t = 5 L/c the advancement of the transmitted probe peak 
is clearly visible: the position of the peak is x pea k ~ 6 L, 
whereas it would have been 5 L if the pulse had prop- 
agated through vacuum (see dashed vertical line in the 
figure) . 

In the unstable regime, for kqL > 7r/2, the instability 
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leads to enormous growth of \£ p (x,t)\ 2 and \£ c (x,t)\ 2 as 
the probe pulse enters the cell. For a spectrally narrow 
pulse centered around 5q = 0, exponential growth sets in 
immediately for that part of the probe which has reached 
the cell boundary at x — 0. For an incoming pulse cen- 
tered around <5o ~ 0.3 c/L the "transient" behavior is re- 
covered. This is illustrated in Fig. [|: at t — the growth 
of the incoming signal has set in in the PCM, but in addi- 
tion there is a clear phase-conjugate pulse with a superlu- 
minally advanced peak traveling to the left. At the same 
time one sees a superluminal "kink" in the transmitted 
probe response: the instability prevents the formation of 
a full gaussian-shaped transmitted probe pulse for this 
set of parameters. The t — lOL/c time plot shows the 
superluminally reflected signal on a larger scale. 
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<(L) <|L) <(L) 

FIG. 5. Unstable (above-threshold) regime of the PCM: 
reflection and transmission of an incoming probe pulse. The 
PCM is located between < x < L (see shaded area). The 
upper (lower) half of each time plot shows the probe (conju- 
gate) pulse |£ p (:r,i)| 2 (- \£ c (x,t)\ 2 ), which is moving to the 
right (left). Parameters used are So = 0.31 c/L, A x = 19.3 L, 
A 5 = 0.144 c/L, kqL = 1.7. 



IV. CHOPPED SIGNALS 

Analytic gaussians are of little value as far as informa- 
tion transfer is concerned, or as a check whether the ob- 
served superluminal effects are in agreement with causal- 
ity |6|,p|,pfj| . For that purpose, one needs an incoming 
modulated or chopped signal. Interesting questions then 
arise: how is the sudden change in the input pulse re- 
flected in the output pulse? How are the reflected re- 
sponse and exponential growth in the unstable regime 
affected by this change? We have already shown else- 
where that for probe transmission in the stable regime, 
the edge of a chopped incoming signal is always transmit- 
ted causally (l5| . Figure || shows the transmitted probe 
response in the unstable regime for an incident gaussian 
pulse which is suddenly switched off at t = 0. 
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FIG. 6. Transmitted probe pulse \£ P [L, t)\ through a PCM 
at x = L for an incident gaussian-shaped chopped pulse (thick 
solid line). Curves a and b correspond to <5o = 0.28 and 
0.3 respectively. The parameters used are the same as in 
Fig. ^| and t tr indicates the time it takes to traverse the cell 
in vacuum. 

We sec that the sudden change in the incoming pulse 
is carried over into the transmitted probe response at 
time ttr = L/c, the time it takes to traverse the cell 
in vacuum. The "information content" of the pulse is 
thus transmitted with the speed of light, after which the 
exponential growth due to the instability immediately 
sets in. Note that the superluminal peak advancement of 
the reflected conjugate pulse remains, just as for the full 
gaussian of Fig. |[ This advancement does not violate 
causality, but is a pulse reshaping effect. The fact that 
the chopped edge of the pulse is transmitted with the 
speed of light can be seen more clearly by expressing 
£ p {x,t) as (App. A) 



£ p (L,t) = J2 

n=0 



dt£ p (Q,t)L n {t,t), (24) 



with 

L n (t,t) = £f (^7 2n _! \k q c^ 
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2„ + :: K C^J(t~t') 2 - {n+ i) 2 r 2 ] ), 
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Here /„ are the modified Bessel functions. £ p (L,t) is 
now expressed as a sum of integrals, in which the n th term 
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corresponds to the contribution after the n round-trip 
time r = 2 L/c. The advantage of this expression is that 
it shows that for an incoming chopped signal, which is 
suddenly switched on at t = to , the transmitted response 
only starts at t = to + L/c, in agreement with causality. 

V. CONCLUSIONS 
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In conclusion, we have theoretically studied the re- 
flection and transmission of wave packets at a phase- 
conjugating mirror. Our main findings are 

(a) In the stable operating regime (for k L < n/2), the 
peak of the transmitted signal can exhibit a superlumi- 
nal peak-traversal time [|l5| . The conditions for this effect 
are that the incoming analytic probe signal be spectrally 
narrow and centered around a frequency sufficiently far 
away from the parametric resonance (5 = in Fig. 0(a)). 
The maximum advancement obtained is ~ 0.88 L/c and 
the maximum ratio (peak advancement) /(pulse width) 
~ 0.08. No such effect is found in the reflected phase- 
conjugate response, whose peak is always delayed wi th re - 
spect to the one of the incoming probe signal (Sec. [II B| ). 
The salient advantage of superluminal peak transmission 
in this stable regime is that the output pulse is undis- 
torted. 

(b) In the unstable regime (for k q L > tt/2), for incident 
probe pulses with a temporal width much larger than 
L/c, a pronounced "transient" reflected phase-conjugate 
response develops before the onset of exponential growth 
due to the instability. This pulse response exhibits a 
superluminally advanced peak, both in phase-conjugate 
reflection and in probe transmission (Sec. II C). The ad- 
vancement and intensity are maximal for kqL close to 
ir/2 and the response is very sensitive to the central fre- 
quency Sq of the incoming pulse. 



(c) We have demonstrated (Sec. IV) that the superlu- 
minal features are observable only for temporally broad 
analytic pulses, in agreement with the principle of causal- 
ity, whereas a sudden (non-analytic) change in the inci- 
dent probe pulse propagates with the speed of light. 

Finally, the question arises how these pulse reshaping 
effects can be observed. For a realistic PCM, consist- 
ing of a cell of length L ~ I0~ 2 m, coupling strengths 
of k c ~ c/L ~ 10 10 s _1 have been reached, so that 
tan 2 («oL) ~ 1 p3]. In order to observe pulse reshaping 
in the stable operating regime of this PCM, one needs 
an incident probe pulse of width ~ 0.1 ns whose peak 
is then transmitted with a superluminal peak transmis- 
sion time ~ 0.01 ns. To enter the unstable regime, the 
PCM has to be operated using pulsed pump beams. The 
pump pulses should be long enough to allow for observa- 
tion of the "transient" pulse response, but short enough 
to avoid the instability effects. Since the width of the 
superluminally reflected and transmitted response is on 
the order of ~ 10 L/c, nanosecond pump-pulse durations 
are required. 



APPENDIX A: TWO-SIDED LAPLACE 
TRANSFORM TECHNIQUE 

In this appendix we briefly outline the two-sided 
Laplace transform (TSLT) technique introduced by 
Fisher et al @. The TSLT is defined as F {r) (x,s) = 
dtF^*'(x,t)e~ st . It is only valid for functions 
F{x, t) that diminish faster than exponentially at times 
t — > — oo jl7]]. Its advantage compared to the the usual 
one-sided Laplace transform is that it also applies to 
functions that do not vanish at t < 0. The starting point 
of the analysis is to apply the TSLT to the four-wave 
mixing equations (||) with F(x,t) — £ p {x,t) or £*{x,t). 
We then obtain the coupled equations 



4=£ p (x,s) + ^£ p (x,s) +ik£ c ^(x,s) = 



dx 
_d_ 



£c,*{ x } s ) ~ f £c,*( x j s ) 4" i^* £p(x, s) — 0. 



(Al) 



Equations ( |A1| ) are solved together with the Laplace 
transforms of the boundary conditions £ p (0, t) = F(0, t), 
where F(0, t) is an incident probe pulse at the entry x = 
of the PCM medium, and £ c (L,t) = 0, so no incoming 
conjugate pulse at the end of the medium. The result is 

£ c (x, s) — h r (x, is) J-"*(0, s) 
£ p (x, s) = h t (x, is) F(0, s), 

with the reflection and transmission amplitudes 

kq sin(/3(L — x)) 



h r (x, S) 



l sin(f3L) + i/3cos(/3L) 



(A2) 



i f3co S (f3(L-x)) + l S in(f3(L-x)) 

ht (x ' = s ■ TlrFVZ a — T^T\ ' (A3) 

|sin(pL) + ipcos(pL) 

where (3 is given by Eq. (^). 

The reflected phase-conjugate pulse £ c (x,t) and trans- 
mitted probe pulse £ p (x, t) at position x in the nonlinear 
medium and time t are then obtained by using the inverse 
Laplace transform. At x = and x = L, respectively, 
they are given by 



£ c (0,t) = ^ 

Z7TZ 



7+200 



dsr C p(is)£p,{Q,s)e st (A4) 



7 — 200 
2 /-7 + 20G 



dst pp (is) fp(0, s) e st 7 (A5) 



7—200 
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where r cp — h r (x = 0) a nd t vv — h t (x = L). The choice M, 
of contour 7 in ( A4 ) and ( |A5[) is in agreement with causal- 
ity, which means in this case that it is to the right of all 
the singularities of r cp (is) (and t pp (is)). The singularities 
in the right-half s-plane give rise to expo nenti al gr owth of 
£ c (x,t) and £ p (x,t). The integrals (A4) and ( A.5) can be 
evaluated by taking the contribution of these poles sep- 
arately and rewriting the remaining integral as a Fourier 
integral, see (0) and (g). and 

They can also be evaluated in another way, which is es- 
pecially insightful when regarding non-analytic, cho pped 
probe pulses (Sec. IV). To that end we rewrite (A.4) and 
(|A|) as 



,.{t,t) = =±f£ (^- 1 / 2n _ 2 [ KQ c^(t-t') 2 -n 2 r 2 
-2BlI 2n \ K() c^/(t-t') 2 -n 2 r 2 



+B™ +1 I 2n+2 KQCy/it-t')*- 



t — t —TIT 
t — t' + flT 



(A14) 



(A15) 



£ c (0,t) = J dt £;(0,t)H c (t,t) 
£p(L,t) = / dt £ p {0,t')H p (t,t), 



(A6) 
(A7) 



with 



1 j-y+ioc , 

Hc(t,t') = — dsr cp {is)e- s{t ^ (A8) 

2m J-y-ioo 

H P (t,t') = — dst pp (is)e- s ^-^ (A9) 

^7Ti J-f — ioo 

In order to evaluate the integral in, e.g., H p (t,t ), t pp is 
first rewritten as the series 



t pp {is) = 2?7^ 



-(n+i)rjT 



(s + if) 



2n+l 



(A10) 



with a = kqc, i] = V s 2 — a 2 and r = 2L/c, the roundtrip 
time of the PCM. One can then easily prove that ( |A10| ) 
is uniformly convergent, which allows for term by term 
integration of flA9|). We use the substitution 



ia ( u A„ 



with 



A n = 



' t-i - {n + \)r 
t-t' + (n + \)t 



1/2 



(All) 



(A12) 



in which t + (n + s)r is the retardation time after (n+ 1) 
round trips. We then arrive at the result Eq. (0). 
Similarly, the conjugate reflected pulse is given by 



£ c (L,t)= I dt S;(0,t)M (t,t) + 

c 

00 r t-nr 

dt£;(0,t')M n (t,t), (A13) 



■2E 



with 



[1] L. Brillouin, Wave Propagation and Group Velocity, 

(Academic Press, New York, 1960). 
[2] C.G.B. Garrett and D.E. McCumber, Phys. Rev. A 1, 

305 (1970); S. Chu and S. Wong, Phys. Rev. Lett. 48, 

738 (1982). 

[3] M. Biittiker and R. Landauer, Phys. Rev. Lett. 49, 1739 
(1982); Th. Martin and R. Landauer, Phys. Rev. A 45, 
2611 (1992). 

[4] A. M. Steinberg, P.G. Kwiat, and R.Y. Chiao, Phys. Rev. 

Lett. 68, 2421 (1992). 
[5] A. Ranfagni, D. Mugnai, and A. Agresti, Phys. Lett. A 

175, 334 (1993). 1089 (1993). 
[6] Y. Japha and G. Kurizki, Phys. Rev. A 53, 586 (1996). 
[7] E. Picholle, C. Montes, C. Leycuras, O. Legrand, and J. 

Botineau, Phys. Rev. Lett. 66, 1454 (1991). 
[8] A. Icsevgi and W.E. Lamb, Phys. Rev. 185, 517 (1969). 
[9] E. L. Bolda, Phys. Rev. A 54, 3514 (1996); E. L. Bolda, 
J.C. Garrison, and R.Y. Chiao, ibid 49, 2938 (1994). 

[10] M. Artoni and R. Loudon, Phys. Rev. A 57, 622 (1998). 

[11] Y. Aharonov, A. Komar, and L. Susskind, Phys. Rev. 
182, 1400 (1969). 

[12] R.Y. Chiao, A.E. Kozhekin, and G. Kurizki, Phys. Rev. 
Lett. 77, 1254 (1996). 

[13] D. Lenstra, in Huygens Principle 1690-1990; Theory and 
Applications, edited by H. Blok, H.A. Ferwerda, and H.K. 
Kuiken (North Holland, Amsterdam, 1990). 

[14] R.A. Fisher, B.R. Suydam, and B.J. Feldman, Phys. Rev 
A 23, 3071 (1981). 

[15] M. Blaauboer, A.E. Kozhekin, A.G. Kofman, G. Kurizki, 
D. Lenstra, and A. Lodder, Opt. Comm. 148, 295 (1998). 

[16] Optical Phase Conjugation, edited by R.A. Fisher (Aca- 
demic Press, New York, 1983). 

[17] G. Doetsch Theorie und Anwendung der Laplace- 
Transformation, (Dover, New York, 1943). 

[18] If the intensities of probe and conjugate waves approach 
those of the two pumps, pump depletion has to be taken 
into account and our analysis is not valid anymore. 

[19] The name "optical quasiparticles" is introduced in view 
of the analogy |l3|pl| with the quasiparticle excitations 
in a superconductor, described by the Bogoliubov-de 



8 



Gennes equations, see P.G. de Gennes, Superconductiv- 
ity of Metals and Alloys, (Benjamin, New York, 1966). 
Just as these superconductor quasiparticles have a mixed 
electron-hole character, the optical quasiparticles consist 
of both a probe and a conjugate part. 

[20] R.Y. Chiao, A.M. Steinberg, and P.G. Kwiat in Quantum 
Interferometry, edited by F. de Martini, G. Denardo, and 
A. Zeilinger (World Scientific, Singapore, 1994). 

[21] H. van Houten en C.W.J. Beenakker, Physica B 175, 187 
(1991). 

[22] M.Y. Lanzerotti, RW. Schirmer, and A.L. Gaeta, Appl. 
Phys. Lett 69, 1199 (1996). 



9 



